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We propose a direct, coherent coupling scheme that can 
create massively entangled states of Bose-Einstein condensed 
atoms. Our idea is based on an effective interaction between 
two atoms from coherent Raman processes through a (two 
atom) molecular intermediate state. We compare our scheme 
with other recent proposals for generation of massive entan- 
glement of Bose condensed atoms. 



Entanglement lies at the heart of the difference be- 
tween the quantum and classical multi-particle world. It 
is a phenomena which facilitates quantum information 
and quantum computing with many qubits. Recently, 
several interesting developments have occurred in studies 
of massively entangled atomic states. Based on the pro- 
posals of Maimer and S0rensen 111, a controlled, entan- 
gled state of 4-ions was successfully created by the NIST 
ion trap group H. Zeilinger and coworkers prepared 
three entangled photon or Greenberger-Hornc-Zcilingcr 
(GHZ) states by selecting from two beams of entangled 
photon pairs || . Entanglement between two atoms and a 
microwave photon was also detected in a "step-by-step" 
process 0]. 

Of these and other related developments, the idea of 
M0lmer and S0rensen [Q is especially interesting. They 
proposed a direct coupling to the multi-particle, entan- 
gled final state through a virtual, intermediate state 
which was a common (quantum) mode of the motion of 
all the ions. Similar type interactions were also proposed 
by Milburn g. Both proposals allow for creation of mas- 
sive entangled states by unitary evolution, starting from 
certain pure initial states. 

A Bose-Einstein condensate of a dilute atomic vapor 
is a convenient source of atoms, well approximated as 
initially being in pure and separable states. S0rensen 
et al. [|| suggested creating massively entangled, spin 
squeezed states from a two component condensate us- 
ing the inherent atom-atom interactions. Spin-exchange 
collisional interactions in a spinor condensate was also 
proposed as a candidate for creating entangled pairs of 
atoms [0J^]. Most of these proposals work in the two 
mode approximation where one motional state is as- 
sumed for each spinor component of condensed atoms. 
Similar studies were previously performed for condensate 
atoms in a double well potential [9|-[T^] . 

Raman transitions mediated by long-range dipole- 



dipole interactions have been proposed in entanglement 
schemes for quantum computing [ p^T^ ]. In this letter, 
we propose a new type of coupling also based on the 
two atom, effective interaction from a Raman process 
through intermediate, molecular states. We show that 
our coupling applied to condensate atoms can achieve 
both massively entanglement, similar to that of M0lmer 
and S0rensen |t[ and improved spin squeezing fel-Mjfl. 
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FIG. 1. (a) The general two-photon Raman scheme via an 
intermediate, excited, molecular transition, (b) Raman cou- 
pling between two spin states illustrated for the atomic and 
molecular case. In the atomic case, single atoms can be pre- 
pared in a coherent superposition of spin up and spin down. 
In the molecular case, pairs of atoms are prepared in a coher- 
ent superposition of spin up and spin down. 

Consider a system (see Fig.[j]) involving two A-type 
atoms whose initial and final states are described by the 
same non-interacting atomic states \g) and \g'). Tran- 
sitions between \g) and \g') for the pair of atoms are 
coupled by two laser fields (frequency u>i, wavevector 
ki,ki — uJi/c, and Rabi frequency f2j) through an in- 
termediate excited state, which is chosen to be a bound, 
molecular excited state. 

We consider an excited molecular state |e&) asymp- 
totically connected, for large internuclear separation, to 
one ground |jtt),/x = g,g' and one excited atom |e). An 
example is the 0~ state, which has been extensively dis- 
cussed in the context of photo-association of ultracold 
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atoms |lq|L7[ ]. The recent experiment by Heinzen and 
coworkers on the production of ground state molecules 
from an atomic condensate by a two-photon Raman 
process |T§| ] provides additional motivation to explore 
our ideas experimentally. The photo-association pro- 
cess 1l6|,|T7lR^ 20 used in Heinzen's experiment relies on 



the transition strength of going from a ground 'free' (two 
atom) to a ground 'bound' (molecular) state via an in- 
termediate, excited 'bound' (molecular) state. On the 
other hand, what we desire is a transition from a ground 
'free' to another ground 'free' state via an intermediate, 
excited 'bound' state. 

The question of whether this is possible or not for 
trapped ultracold cloud of atoms does not seem to de- 
pend on the sample density (in the weakly interacting 
limit), but will depend, to a large degree, on the trap 
strength and the excited, bound state structure. That 
is, we can determine the strength of the transition by 
considering the collision of a pair of atoms in the com- 
bined trap and molecular potentials and then summing 
over the number of pairs of atoms available. 

Assuming all trapping potentials for the ground and 
excited states to be harmonic, there is a separation of 
the trapping potentials in terms of the center-of-mass R 
and relative coordinates f of the two atoms. When short 
range interactions are approximated by their optical con- 
tact forms, Wilkens and coworkers provide analytic so- 
lutions for the center-of-mass and relative motion of two 
interacting atoms inside a harmonic trap [ pl| . In prin- 
ciple we can find all bound states in the r potential due 
to the external trapping potential, even with a realis- 
tic molecular interaction potential [p2| . In general one 
can include all relative motional states and still be able 
to find the effective coupling between the selected elec- 
tronic states p3| . We consider, however, only the lowest 
unbound state of the relative motion |0) in the electronic 
ground state, a situation well approximated by Bose con- 
densed atoms. 

If we choose a state fairly deep in the excited state 
molecular potential, then adjacent molecular states will 
be well separated in energy and we can consider coupling 
to only one intermediate bound state with vibrational 
quantum number mj. Furthermore, for a sufficiently 
large detuning of the coupling lasers from the excited 
state eb, {mi,), we can perform an adiabatic elimination of 
the excited state. Neglecting configurations not directly 
involved in the two photon process p4[ , the Hamiltonian 
for the two atoms initially in the g state can be written 
as 
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V tR (R) + 2foj Mfl + E 



(1) 



Vr = ^Imrnfe^e^-x^KOfag-A + h.c. (2) 



Auj = u>2 — u>i, K = ki — &2, VtR is the trapping po- 
tential for the center-of-mass motion, hoj^ g is the energy 
difference between the atomic states g and g or g' and, 
for simplicity, we have assumed the same relative mo- 
tional state |0), with energy E , for the pair of atoms in 
g or g' (i.e. atoms in state g or g' see the same trapping 
potential). Since |e&) is asymptotically connected with 
\e)\fi) + \n)\e) (for /i = <?, g') we can set the dipole matrix 
element di = d! i = d 7 for the to |e) transition of atom 
i. Then Q,r = fl 1^3/2 A is simply the two-photon Rabi 
frequency coupling between atomic states g and g' . Typi- 
cally \ki\ ss |fc2| and we can use r?om b = (0| cos(/ci • r)|m{,) 
for both free-bound transition amplitudes. |??om b | 2 is es- 
sentially the Franck-Condon factor, or a measure of the 
strength of the free-bound transition and will vary con- 
siderably with r for small values of r, due to the relatively 
short range of the molecular interactions. It is therefore 
important to pick the intermediate state |m&) which re- 
sults in a large l^omj- These values can be determined 
for selected molecular states from the results of photo- 
association experiments. They can also be computed di- 
rectly if accurate potentials are available. Detailed dis- 
cussions are given in 1 17 ll||2(| . 

For two counter-propagating waves (K /2 ps k\ ~ — k%) 
and fj, — g, the above result is similar to the single 
atom case, which is often referred to as Bragg diffrac- 
tion [psj . The coupling involves the simultaneous ab- 
sorption and stimulated emission of photons, i.e. the 
elementary process involves a total of two photons. In 
contrast to the single atom case, a pair of atoms are now 
involved. Hence for Bose condensed atoms with w 0, 
Bragg diffraction produces pairs of atoms in the state 
a\pi ps 0,p 2 ~ 0) + /3\pi = HK/2,p 2 = UK/ 2). The mo- 
mentum shift per atom is only half the value of atomic 
Bragg diffraction, due to the fact that only two photons 
are involved for the pair of atoms. Hence the resonance 
(energy conservation) condition Aui occurs at half the 
atomic Bragg resonance. 

For a Raman process with two co-propagating waves 
(ki ps k2, K/2 ps 0) between two nearly degenerate states 
g and g' (as depicted in Fig.^(b)), the resonance con- 



dition is Auj 



2oj 



twice the atomic resonance; i.e. 



where 



the Raman process changes each atom's state (g — ► g') 
and so each atom acquires the energy deficit Tiuj g r g . Note 
that, in this case, there is no R dependence in Vr. 

We now consider many atoms in a spin 1/2 system to 
investigate the extent to which Raman coupling via an 
intermediate molecular state can produce spin squeez- 
ing [^5) and correspondingly massive entanglement ||. 
If we designate \g) and \g') as | T) and | J.) for spin up 
and down, respectively, our Raman coupling for the two 
particle case is of the form 
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[(ITXil)i®(ITXil)2 + (UXTI)i®(UXT|)s 



with GHZ-type state |GHZ)jv at selected times (see Fig. 



(3) 



where u„ =a . „ z are the Pauli matrices. This is a new 
type of coupling not widely discussed before. Its ability, 
however, to generate entanglement should be obvious as 
two atoms perform conditional evolution at all times. 

We first distinguish our coupling scheme from other 
relevant models. In the original scheme of M0lmer and 
S0rensen M, the two atom coupling takes the form 
ifiQiia^ <8> <Jx^ , a form different from ours. If the 
pair-wise interaction acts indiscriminately for all pairs 
of atoms, it is convenient to analyze the effect of such 
couplings for many atoms in terms of a collective spin 
operator = X^ "^ jA* = x,y,z, where the sum is 
over the number of atoms N |ioUl2|,j26 27 1. One can 



then show that the Hamiltonian for the M0lmer and 
S0rensen scheme becomes Vm = J2i<j \ fi£lR&x ®Vx — 
jhQff,{Jx — N). In contrast, our scheme gives 
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(J'-Jy)- (4) 



Recently, a many body, two mode coupling scheme was 
proposed by S0rensen et al. M. They considered a two 
component (i.e. | |) and | |)) condensate weakly inter- 
acting via s-wave collisions described by a mean-field. In 
the approximation where each component has the same 
spatial mode, the interaction in terms of the collective 
spin operators J x ,y,z is of the form Vs = ^hrijiJ 2 . 

In order to compare the various coupling schemes, 
we investigate the time evolution assuming a pure ini- 
tial state with a fixed, total number of atoms N. Us- 
ing the notation of second quantization, where and 
a,j are the operators for creating and annihilating a 
particle in state j (j =t,l), the collective spin oper- 
ators can be written as J x = (aja-f + ajajJ/2, J y = 
«(ojof — ajajJ/2 and J z = (aja-f — atajJ/2. To nu- 
merically calculate the time evolution, we expand the 
wave function as = X)m=o c m(t)\m) jJiV — m)t, where 
\m)j = (aj) m |vac) / \fm\, with the initial conditions given 
by the c m (0)s. 

For N an even integer ]2S| ], the time evolution opera- 
tor for Vm takes a simple, analytical form at ^Ipt = tt/2, 
producing a massive GHZ-type wave function | GHZ) at = 
[|iV)||0)x + ?7|0)-f|iV)jJ/\/2, where rj is purely a phase 
factor. In the context of Bose-Einstein condensation, 
this is an example of an interesting fragmented conden- 
sate While the M0lmer and S0rensen coupling (Vm) 
produces perfect GHZ-type states at selected times, the 
S0rensen et al. spin squeezing scheme (Vs) and our 
scheme (Eq.([|)), in general, do not produce exact GHZ- 
type states. From numerical simulations, however, we 
find that our model can produce more than 50% overlap 




FIG. 2. The time dependent probabilities |co| 2 (solid line) 
and |cjv| 2 (dashed line) of being in states |JV)jj0)j and 
|0)||iV)j, respectively, for (a) the M0lmer and S0renson cou- 
pling (Vm) and (b) our "molecular Raman" coupling (Eq^). 
The initial conditions are c m (0) = 8o m . The number of atoms 
in the simulation is N = 10 3 . The projection onto the [GHZ) n 
state corresponds to | Co | 2 = jcjv| 2 . In (a) it is apparent 
that the coupling Vm produces a perfect | GHZ) at state at 
Q.Rt — 4-7T. While in (b) it appears that at Slut pa 0.28 (in- 
dicted by the arrow) our coupling scheme results in about 
50% of the atoms in the |GHZ)jv state. 

We can also compare the achievable spin squeez- 
ing between our scheme (Eq.(^)) and that of S0rensen 
et al. (Vs), using the squeezing parameter £ 2 = 
N(AJ Hl ) 2 /((J n2 ) 2 + (Jff 3 ) 2 ), where n u i = 1,2,3 are 
mutually orthogonal unit vectors [pp5|] . 

S0rensen et al. |(| have shown that £ 2 (i > 0) < 1 for 
some set of n,'s. Their scheme (Vs) is in fact the one- 
axis twisting model considered by Kitagawa and Ueda 
earlier [H . In this case the problem can be solved ana- 
lytically, with the result (AJiJ min ~ 7V 1/3 . 

On the other hand, our coupling resembles the the two- 
axis countertwisting model of Kitagawa and Ueda |l5| , 
and has to be solved numerically. In the limit of large 
N and with the condensate initially in one spin state, 
J z = -N/2, one can show that (AjiJ^ sa 1/2. The 
optimal squeezing in this case occurs along x + y. This 
result can be easily verified by making a semi-classical 
approximation in the dynamical equations for J x and 
J y . We find that the time scale of reaching maximum 
squeezing is ~ 1/(NQ,r) (see also Q). For condensates 
containing 10 6 atoms, even with a very weak coupling 
fl R sa 1 (Hz), the maximum squeezing is reached within 
a microsecond. 

In Figure |3|, we show the numerically computed, min- 
imum squeezing parameters, £ 2 , as a function of time 
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for the coupling from our scheme (Eq.(|j)) and from the 
S0rensen et al. scheme (Vs). In contrast to the S0rensen 
et al. scheme, our scheme achieves better squeezing at 
an earlier time. In addition for Vs, the direction ni 
along which minimum squeezing is observed varies with 
time H[l5| . While for the coupling of Eq. (Q) , it is fixed 
along x + y jl5| . 




FIG. 3. The time evolution of the minimum spin squeezing 
parameter £ 2 for (a) the coupling Vs and (b) the coupling 
given by Eq.(§). 

In order to realize our coupling scheme experimen- 
tally, the molecular coupling (Fig. 1(b), right-side) needs 
to dominate over the atomic coupling (Fig. 1(b), left- 
side). That is, we need to achieve a molecular coupling 
Q% = OiOH^lVAjwf » £1% = nifi5/A Aj where 
Am /a ar e detunings from the molecular and atomic in- 
termediate states, respectively. We should also have 
Am^Tm, the excited molecular state linewidth, to min- 
imize spontaneous emission which would lead to decoher- 
ence and loss of atoms. These constraints suggest that 
a deep molecular bound state with significant transition 
strength should be chosen in order to maximize the de- 
tuning from the atomic transition and achieve a sufficient 
detuning from excited molecular states. In addition, if 
Aa 3> Am further suppression of the atomic transition 
might be possible by a suitable choice of laser polariza- 
tions. 

Another mechanism for suppressing the atomic cou- 
pling with respect to the molecular coupling is the two- 
photon, Raman resonance detuning. For Bragg diffrac- 
tion, the transition via a molecular coupling occurs at 
half the detuning of the transition via an atomic coupling 
(50 kHz vs. 100 kHz, respectively for sodium). If the Ra- 
man transition involves changes in the internal state of 
the atom, then the frequency for the Raman transition 
via a molecular coupling will be at twice the frequency 
of the Raman transition via an atomic coupling (see Fig. 
1(b)). For Raman transitions between Zeeman sublevels 
in modest magnetic fields this frequency difference can 
be several MHz, which would greatly suppress the single 
atom transition. 

There are several advantages of our coupling scheme. 
The different Raman resonance frequency is a clear signa- 
ture for the Raman transition via a molecular intermedi- 
ate state. In addition, for Bragg diffraction via the molec- 
ular coupling, the atoms would move at half the speed of 



atoms that have undergone atomic Bragg diffraction. 

Our scheme is based on an engineered interaction that 
can be turned on and off like the scheme of M0lmer and 
S0rensen (Vm). In contrast to the S0rensen et al. cou- 
pling scheme (Vs), our scheme works for non-interacting 
ground state atoms, which can decrease the noise due 
to atom-atom interactions in a U(l) symmetry break- 
ing condensate state. Our scheme also achieves the same 
level of squeezing and the same high value of overlap with 
the massive GHZ state, even as the number of atoms is 
increased. 

Finally, we note that our discussions above can also 
be applied to a pair of different species of atoms. For 
example, Raman coupling using a molecular intermediate 
state of the Li-Cs dimer could create entangled pairs of 
Li and Cs. 
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